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Forward dihadron correlations in the Gaussian approximation of JIMWLK 



We compute forward dihadron azimuthal correlations in deuteron-gold collisions using a Gaussian approximation 
for the quadrupole operator. The double parton scattering contribution is found to be part of our dihadron calculation. 
We obtain a good description of the PHENIX data for the azimuthal-angle dependent away side peak and a relatively 
good estimate for the pedestal contribution. 
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1. Introduction 

At high energy or, equivalently, small x, the interactions of hadrons are expected to be dominated by nonlinear 
strong color fields. A convenient effective theory approach to studying these color fields is provided by the Color 
Glass Condensate (for a review see e.g. ||T|). 

Measurements of correlations between two forward hadrons in dAu collisions at RHIC |'2','3 1 seem to show indica- 
tions of "initial state" or "cold nuclear matter" effects that are significantly stronger than in pp collisions or at central 
rapidities. The upcoming LHC proton-lead collisions will provide more opportunities to study these phenomena in a 
wider kinematical range. 

These observations have provided an impetus for renewed interest in the gluonic correlations included in the 
JIMWLK evolution, see e.g. [4]. In particular it was argued that the result of a full JIMWLK evolution, at finite Nc, 
can quite accurately be captured by the so called Gaussian approximation, relating higher point Wilson line correlators 
to the two-point function. These recent theoretical developments were not fully reflected in the pioneering calculations 
of dihadron correlations in |5 6 1 where a factorized approximation is derived in a certain kinematical limit. The main 
purpose of this work (see Ref. |7|) is to implement the Gaussian approximation, which so far has only been tested for 
particular coordinate space configurations, in a full calculation of the dihadron correlation. 

2. Single inclusive baseline 
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We shall here use the dipole operator obtained from solving numerically the BK evolution equation using the 
Balitsky running coupling prescription (SI. The BK equation requkes an initial condition, for which we take the 
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Figure 1 : n production in forward p+p collisions com- 
pared with BRAHMS |12 | and STAR |[13] n° data 
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Figure 2: Nuclear suppression factor 7?dAu compared 
with PHENIX 13 and STAR 1 13 1 tt" data. Note that the 
STAR data is at a slightly different rapidity. 
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with 2j = 0.2 GeV^ for the proton at xq = 0.007 similarly as in Ref. fTO\ . Following Ref 111] we can then calculate 
single inclusive hadron production in forward rapidities. The comparison with STAR 7r° and BRAHMS charged 
hadron data is shown in Fig.[T] This parametrization overestimates the pion yield by a factor ~ 3, but describes the 
charged hadron yield correctly. We prefer to use the same parameters as Ref. 1 10 1 to see clearly the effects of including 
the Gaussian quadrupole. We expect most of the overall normalization error to cancel in the per trigger correlation 
and not affect the systematics of the away side peak. 

For a nuclear target we use larger initial saturation scale, and we set Q^^ to 0.72 GeV^ at same xq = 0.007 
which fits the most central PHENIX R^au data shown in Fig.|2] However it seems quite impossible to simultaneously 
describe STAR minimum bias data. We consider the normalization issue of the single inclusive baseline as the largest 
uncertainty in our calculation. 

3. Dihadron correlations 

Consider a large x quark with momentum from the probe deuteron or proton, propagating eikonally through 
the target nucleus or proton. It can radiate a gluon with momentum - zp^ and is left with longitudinal momentum 
q'^ - (1 - z)p^- In the high energy limit the scattering of both the quark and the gluon can be described by the 
eikonal approximation, where they pick up a phase given by a Wilson line in the color field of the target. The detailed 
derivation of the double inclusive cross section is performed in Ref. H) and results in the following expression for the 
qA — > qgX cross section: 

dor^^^"'' _ . . r d^x, dV, d^b, d2bJ-^,,^.(^,_^^,^,^^.(,,_,^, 



,.C,6(p -k -,)J — — — — 



dk+ (fkr d^+ d^qr " " ^ ' j (2n)^ (In)^ (Itt)^ (liif 

X "^^M^'r " ^'T)^U^-r - br) {5<''(br,xr,b'r,x^) - 5<^'(br,xr,z'3-) - 5<2>(zr,x^,b^) + ^(^^(zr, z'^-)} , 

with Zj- = zx/- + (1 - z)b7- and likewise, ilj - zx^ + (1 - z)h'j. The momenta of the produced gluon and quark are 
kj- and qj- respectively. Likewise, X7^,x' should be interpreted as the transverse position of the gluon and of 
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Figure 3: The quark-gluon parton level correlations 
near forward RHIC kinematics Shown are the "naive 
large-A^c" approximation (labeled as "elastic" in the plot) 
used in [6| and our Gaussian approximation for the 
quadrupole and its large-A^j, limit. Note that fixed A(p- 
independent pedestal of 0.0002 GeV '* is added to the 
"elastic" approximation for the purposes of visualiza- 
tion. 



Figure 4: The n azimuthal angle correlation compared 
to the PHENIX |3| dAu result. The Acp independent 
pedestal is adjusted to fit the experimental data. The ini- 
tial saturation scale is gso - O-'^^ GeV^ at xq - 0.007 
and the large-A^c limit of the Gaussian approximation is 
used. 



the quark, in the amplitude and the complex conjugate respectively. Here 5'"' are operators constructed from Wilson 
lines whose expressions can be found e.g. from Ref. fS], and is the q ^ qg splitting function. We set the scale at 
which the parton distribution functions, fragmentation functions and coupling constant are evaluated, to be equal to 
the transverse momenta of the trigger hadron. 

In Fig. |3]we show the parton level dihadron production cross section obtained using different approximations 
for the quadrupole when evaluating S^'^K We notice that the width and especially the height of the away side peak 
are modified when replacing the naive approximation, used e.g. in |6|, by the Gaussian approximation [14J. The 
comparison with experimental data is shown in Fig. |4] where we compute the coincidence probability measured by 
the PHENIX collaboration (3]. Once we have adjusted the A<p independent pedestal contribution we obtain a good 
description of the A<p dependent part of the data. 



4. Double parton scattering 

The Aip independent pedestal in the experimentally measured coincidence probability is mostly due to the double 
parton scattering contribution where two hadrons are produced independently of each other. In previous works ll6l [T5l 
this process has been considered to be completely separated from the correlated dihadron production. However, when 
one uses the Gaussian approximation for the quadrupole function one can show that Eq. (|2| actually contains part of 
this DPS (Double Parton Scattering) contribution. 

The DPS contribution is obtained when gluon is emitted far away from the quark: Ibj - b^| ~ \xt - x'j.| ~ l/Qs, 
uj = Ibj- - X7I » l/2s- In this limit the only surviving correlation in 5^"^* is the product of an adjoint representation 
dipole at the location of the gluon and a fundamental representation one at the location of the quark. We call this 
kinematical regime the "DPS" limit. In this limit the Wilson lines of the quark and the gluon are uncorrected and we 
have, using the fact that the expectation values must be color singlets, 

5<4)(br,xr,b;-,x'r) ~ 5<,'^s(br,xr,b;^,x^) = (D(br, b^)) /D2(xr,x;^) - -^V (3) 
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where D is a correlator of two Wilson lines (two point function) in the fundamental representation, and the expectation 
value of - 1 /N^^ can be identified as the two point function in the adjoint representation. 

For massless quarks the four-point function 5'*"^' causes the integral ^ to diverge logarithmically in the DPS 
limit. Physically this means that the quark emits a very small transverse momentum gluon. The quark and gluon 
subsequently scatter independently off the target. The DPS limit corresponds to a splitting happening a long time 
before the interaction with the target. This logarithmically divergent contribution must be regulated by confinement 
scale physics in the wavefunction of the projectile. It is in fact exactly the kind of contribution that is represented by 
double parton scattering, discussed in this context e.g. in Ref. (16]. 

We add to the dihadron production cross section the DPS contribution. In order to avoid double counting the 
logarithmically divergent DPS part is subtracted from Eq. ([2]). This forces us to introduce an arbitrary (soft) cutoff 
scale ~ ^QCD- contribution is then divided into two separate parts. The first one corresponds to taking two 

partons from the same nucleon in the deuteron, described by a single-nucleon double parton distribution function, 
which we model by implementing a kinematical constraint x, + xj < 1 following Ref. lfT6l . The second contribution 
involves taking one parton from the neutron and the other one from the proton, which is not bound by the same 
kinematical constraint. Once we have obtained the corresponding parton distribution functions, the DPS yield is 
essentially the single inclusive yield squared. 

As a result we obtain an order-of-magnitude estimate for the Acp independent pedestal. When comparing with 
PHENIX data, we obtain for the trigger transverse momentum range 1.1 ... 1.6 GeV a pedestal 0.11 GeV whereas 
the experimental value is 0.176 GeV Similarly for the trigger transverse momentum 1.6... 2 GeV we obtain 
0.08 GeV \ and the experimental value reads 0.163 GeV ' ISj. 

5. Conclusions 

We have shown that the previously used "naive large-A^c" approximation for the quadrupole is not accurate, and 
the width and especially the height of the back-to-back peak is modified when more accurate Gaussian approximation 
is used. In addition, the naive approximation misses an important logarithmically divergent DPS contribution, which 
is included consistently in our work. 

We obtain a relatively good estimate for the A(f> independent pedestal background, and a good description of the A(p 
dependent part of the PHENIX data. We also point out that the "naive large-A^c" approximation clearly underestimates 
the height of the away side peak. 

We thank K. J. Eskola, I. Helenius, R. Paatelainen, B. Schenke, M. Strikman and R. Venugopalan for discussions 
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of Particle and Nuclear Physics. This work has been supported by the Academy of Finland, projects 141555 and 
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